Cooper pair current in the presence of flux noise 
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We study the efTect of the flux noise on the Cooper pair current and pumping. We generahze the 
definition of the current in order to take into account the contribution induced by the environment. 
It turns out that this dissipative current vanishes for charge noise but it is finite in general for noise 
operators that do not commute with the charge operator. We discuss in a generic framework the 
effect of flux noise and present a way to engineer it by coupling the system to an additional external 
circuit. We calculate numerically the pumped charge through the device by solving the master 
equation for the reduced density matrix of the system and show how it depends on the coupling to 
the artiflcial environment. 



INTRODUCTION 



state during the whole evolution .123111 



Geometric phases are ubiquitous in quantum physics. 
Both the Abelian Berry phasc"'^ and the non- AbeliaiP 
generalization were theoretically discovered in the eight- 
ies. During the last decade, there has been a renewed 
interest for their possible use as tools to manipulate quan- 
tum information. The so-called geometricP or holonomic 
quantum computatiorP relies on the fact that it is possi- 
ble to build unitary transformations which depend only 
on geometric properties of an abstract parameter space. 
The main advantage of this approach is that such geo- 
metric tran sformations are robust against certain types 
of noise.EH 

While the Berry phase has been observed in many dif- 
ferent quantum systemsP^^ clear experimental evidence 
of the non-Abelian adiabatic geometric phases is still 
missing despite different proposals for their theoretical 
implementation.ESli^ One of the main reasons for this 
is that the system must be cyclically steered in an adi- 
abatic fashion. In other words, the time to implement 
a geometric non-Abelian transformation is long with re- 
spect to the dynamical time scale of the system. This 
renders it difficult to protect the system from environ- 
mental noise and thus decoherence can become an issue 
in the implementation of geometric quantum computa- 
tion. 

A step towards understanding the effect of the environ- 
ment and the robustness of steered quantum evolution 
has been taken in Refs. [TTHTO] where it has been shown 
that ground-state evolution is robust against relaxation 
and dephasing induced by a low temperature dissipative 
element. 

Among different physical systems, the Cooper pair 
sluic^2Sl ig a, promising candidate for studying the en- 
vironmental effects on geometric phases. In fact, the 
pumped charge through the Cooper pair sluice is known 
to be directly connected to the geometric phases .l^^^^^ll 
This has allowed the measurement of the Berry phase in 
such a system.[2l This has paved the way for ground-state 
geometric quantum computating in which the quantum 
information is manipulated through a geometric opera- 
tor but the system is kept in a doubly degenerate ground 



Furthermore, transmon quibits coupled to a supercon- 
ducting cavity offer potential systems for the observation 
of the non-Abelian phases-l^fil They constitute the other 
superconducting platform in which the Berry phase has 
been observed experiment ally.'^'^ZIlll 



The Cooper pair sluice operates in the charging regime, 
i.e., the charging energy dominates over the Josephson 
energy. For this reason, the noise induced by charge fluc- 
tuations is likely to be dominant and it has been already 
studied in Refs [T7| and [T51 In this paper, we analyze 
the effect of another kind of noise: the flux noise. This 
can be induced, for example, by a fluctuating magnetic 
field. We find that the definition of the current must be 
extended to include the contribution induced by the en- 
vironment. This dissipative current is directly related to 
the symmetry of the system and can be calculated using 
the master equation approach. 

As a practical application, we discuss the effect of flux 
noise on Cooper pair pumping. Instead of a purely the- 
oretical approach, we study a possible implementation 
in a realistic experimental setup. We couple the Cooper 
pair sluice to an artificial noise source produced by an 
external circuit. The main advantage of this approach 
is the possibility to control in situ the coupling strength 
between the system and the noise source. 

The article is organized as follows. In Sec. |lTj we derive 
the general expression for the current operator, starting 
with the emphasis on the symmetries of the problem and 
then on the connection with the master equation. In Sec. 
|TTTj we apply these results to the pumping process, and 
identify the different contributions to the current and 
the pumped charge due to the environment. We analyze 
the particular case of flux noise which perturbs the phase 



across the superconducting loop. Section IV presents the 
circuits to engineer the artificial fiux noise environment. 
In SecjVj we discuss how the pumped charge is influenced 
by the artificial noise. 
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II. DEFINITION OF THE CURRENT 
OPERATOR 

Let us assume that a quantum system is in a state 
described by the density matrix p. The expectation value 
of a quantum observable A is 

(i)=Tr(M)-E(*l^^l*)' W 

i 

where the last term is the explicit expression in a time- 
independent basis If A is time-independent, the 
time derivative of the expectation value is 

i ^ ^ 

Using the von Neumann equation we obtain the well- 
known Ehrenfest theorem 

|(i>=-^Tr(p[i,ff]), (3) 

where H is the Hamiltonian which determines the dy- 
namics. Having in mind the usual relation between the 
charge and the current, we can say that the above equa- 
tion defines the current associated with the A operator. 

The obtained result is valid for a general quantum sys- 
tem. Let us consider the case in which the total sys- 
tem is composed by a subsystem S and the environ- 
ment E. We are interested in the dynamics and ob- 
servables related to the subsystem S. To simplify the 
discussion, in the following, we will refer to S simply as 
the system. The total Hamiltonian can be written as 
H = Hs + He + Hi, where Hs denotes the Hamiltonian 
of the system. He is the Hamiltonian of the environment 
and Hi describes the interaction between the system and 
the environment. The A operator acts only on the system 
space, i.e., [He, A] = 0, and, from Eq. (jsj) we obtain 

j^{A)^-^{Tiip[A,Hs]) + Trip[A,Hi])}. (4) 

The trace in Eq. Q is over all the degrees of freedom and 
can be split into the trace over the degrees of freedom of 
the system and the environment: Tr(y4) = Tr5[Tr£;(A)]. 
By noticing that both A and Hs act only on the system 
degrees of freedom and that TiEip) — ps (where ps is 
the reduced density operator of the system), we obtain 

f^{A)^-^{Trsips[A,Hs]) + TTip[A,Him. (5) 

The two contributions in the current associated with A 
have different origins. The first term on the right of Eq. 
^ resembles the current for an isolated system. In fact, 
if the system does not interact with the environment. 
Hi = 0, this is the only contribution present. The second 
term can be interpreted as an additional contribution to 
the current due to the interaction with the environment. 



In the following, we refer to the two contributions as the 
system current and the dissipative current, respectively. 

The case in which Hi commutes with the A operator 
is particularly interesting. Here the dissipative contribu- 
tion to the current vanishes. However, the interaction 
with the environment influences the evolution of ps and 
thus it can modify the current. 

A. Connection to the master equation 

Equation ^ is written in order to emphasize the sym- 
metry of the problem. However, if the system and the 
environment have no particular symmetry, it cannot be 
used in general to estimate analytically the contribution 
of the dissipative current since the calculation involves 
the full density operator. 

The dynamics of the reduced density matrix of the 
system ps is obtained by writing the formal solution of 
the von Neumann equation, expanding it in powers of the 
system-environment coupling and then taking the trace 
with respect to the environmental degrees of frecdom.l^Sl 
With this procedure we arrive at a standard form of the 
master equation 

where the dissipative term is denoted by C. 

We can apply this formalism to the case discussed 
above. As assumed in the derivation of the mas- 
ter equation, Tte{^) = ^ and, thus, Tr(^i) = 
Tis{TrE{^A)} = Trs(^i). Using this result, Eq. 
and Eq. @, we have 

= -jTrs{[Hs,Ps]A) + Trs{CA). (7) 

The first contribution on the right matches with the sys- 
tem current in Eq. ([s]) and, hence, we have an explicit 
expression for the dissipative current 

-'-TTip[A,Hi])^TTsicA). (8) 

This equation allows us to calculate the dissipative con- 
tribution to the current using the expression of C from 
the master equation for the reduced system dynamics. 

Equation ([s]) has important implications involving the 
quantities that are preserved during the evolution in pres- 
ence of environmental noise as recently discussed in Ref. 

1301 



III. APPLICATION TO CHARGE PUMPING 

The above analysis helps us to set the framework to 
discuss the pumping process in the prese nce of an envi- 
ronment. We consider a Cooper pair sluicd^SElini! shown 
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in Fig. [T] It consists of a superconducting island sepa- 
rated by two superconducting quantum interference de- 
vices (SQUIDs) with controllable effective Josephson en- 
ergies Jl,r- The electrostatic potential on the island can 
be controlled by varying a gate voltage Vg. The experi- 
mental access to the parameters Jl,r. and Vg allows for 
a full control of the quantum system and makes it an ex- 
cellent prototype for different applications. Several steps 
have been taken with the study of the connection be- 
tween Cooper-pair pumping and geometric phases, both 
in its AbelianSElJ and non-Abelian versi orP^" ^, the ro- 
bustness of the ground state pumpin^i^tiSlllHlIl and ge- 
ometric Landau-Zener-Stiickelberg interferometry.'^ 

The Hamiltonian of the sluice Hs is the sum of the 
charging Hamiltonian Hch ~ Ec{n^ngY and the Joseph- 
son Hamiltonian^^ ■^^ 



-JlCos [ - 



Jr cos ( - 



(9) 



where (p = ipn + '.f'L is the superconducting phase dif- 
ference between the two leads, Ug = CgVg/{2e) is the 
normalized gate charge, Eq = 2e^/Cs is the charging 
energy of the sluice, Cg is the gate capacitance, and 
Cs the total capacitance of the island. We denote with 
fik — —idip^. (k = L, R) the Cooper pair number opera- 
tor of the fcth SQUID, and with 9 = {<p>R — (pl)I1 and 
n = —ide the operators for the superconducting phase 
and the number operator of excess Cooper pairs on the 
island. If the device operates in the charge regime, i.e., 
Ec ^ max{ Ji,Jfl} — Jf^'p and the gate parameter is 
close to half integer, only the two lowest energy charge 
states are important for the dynamics and we can adopt 
the two-state approximation. Let |0) and |1) denote the 
states with no and one excess Cooper pair on the island, 
respectively. For the system in consideration, it is conve- 
nient to reduce the Hilbert space restricting it to states 
with well-defined ip. In this case, p = 27r$/$o can be 
treated as a real number and it is determined by the 
magnetic flux through the large superconducting loop 
in Fig. pa). 

Starting from the definition of the charge operator 
through the fcth SQUID Qk — —2enk, the discussion in 
Sec. [n] can be interpreted in terms of physical quanti- 
ties. If the dynamics is influenced by the environment, 
the average current through the kth SQUID is 



The second contributions in Eqs. ( 10 1 and ( 11 ) represent 
an additional dissipative current 



lip - - 

{Itn = -YTTip[nk,Hi]) = TvsiC Qk) 

induced by the environment. 



(12) 



It is convenient to write Eq. (11) in the eigenbasis of 
Hs- Let \g) and |e) be the eigenstates oi Hs and M„j„ = 
{m\M\n) with m,n = g,e for any operator M. Equation 
( 11 1 becomes 



{-^k} PS^gg-^k,gg PS,ee-^k,ee ~^ {Qk.gg Qk,ee)^gg 

+ 2^e{ps,geh.eg) +2^e{Cg,Qk,eg), (13) 



where we have used the fact that, from the symmetries 
of the master equation, Cee — ^^gg and Ceg — C*g^. We 
identify the first line as the dynamic current and the sec- 
ond one as the geometric current. Both are composed 

of standard contributions ~ Ps,gglk,gg + PS.eeIk.ee 

and = 2^e{ps,gelk,eg), and dissipative contributions 



Z),diss 



{Qk,gg Qk,ee)^gg and Z^;^ — 2die(^/lgf^Q f^^^g j 

The pumped charge through the fcth SQUID is defined 



as 



= 2 



'^e{ps,gelk,eg)dt + 2 {CgeQk,eg)dt, 



(14) 

and it is also composed of a standard and a dissipative 
part. From the above expressions we can write the cur- 
rent operator for the average current across the device 
/ = {II + Ir)/2 and the corresponding average geomet- 
ric current in the sluice. 

Note that the dissipative contributions jj^-'^^^^ and 
jG',di!,s typically small with respect to ijp and I^. 
If A is the effective coupling between the system and the 
environment, the dissipative contribution in the master 
equation £ scales as A^. Since the master equation ^ is 
derived in the limit of weak coupling between the system 
and the environment, we expect a small contribution in 



Eq. ( 13 1 from the dissipative currents. However, there 



can be cases in which, the dissipative contributions could 
be detectable if we reduce and I^. 



A. Charge noise environment 



j^{Qk)^j-Tys{ps[nk,Hs]) 



2ie 



Ti{p[nk,Hi]). 



(10) 



The current operator for a closed system is usually de- 
fined as Ik — [fik , Hs] and it corresponds to the first 
term on the right side in Eq. (10) (see Ref. 22 ). Using 
Eq. (Isl), we have 



(/fe) =Tr5(/5s/fe)+Trs(£ Qu). 



(11) 



When the system is in the charge regime, the main 
source of noise originates from the fluctuations of the 
gate voltage. References P^THTO] focus on the study of 
this charge noise induced by the environment. If we 
write the Cooper pair number operator and the charge 
on the island as n = — fiji and Q = —2en, re- 
spectively, the corresponding system-environment inter- 
action is Hi = —2eXn (g) 5V, where A is the system- 
environment coupling const ant an d SV acts on the envi- 
ronment degrees of freedom.lSMl The operators fik com- 
mute with h and, thus, there is no dissipative current 
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FIG. 1: (Color online) (a) Circuit diagram of the Cooper pair sluice (on the left) and the artificial environment circuit (on 
the right). The experimentally controlled parameters are Jl,r and Ug, and = 2n^ / Phio is the phase difference across the 
device and "l>o is the flux quantum. The system is coupled to the artificial environment through mutual inductance M. The 
engineered environment is composed of a resistor R associated to a voltage source SV, an inductor L and an experimentally 
controllable SQUID. The SQUID allows us to change the effective noise spectrum of the environment, (b) Time dependence of 
the parameters Jl,r and Ug during a Cooper pair pumping cycle. 



i^k^^^) through the fcth SQUID and no total dissipative 
current contribution. As discussed in Sec [llj even if the 
dissipative current vanishes, the total current is modi- 
fied by the interaction with the environment through the 
dynamics of ps [see Eq. (10)]. 

Since no dissipative current can be induced by charge 
noise, i.e., [Q,Hi] = 0, Eq. (jsj) requires that the term 
Tr5(>C Q) calculated from the master equation vanishes. 
This can be verified by a direct calculation using the 
master equations as discussed in Appendix |X] Hovifever, 
this result depends critically on the form of the master 
equation used and the approximations done. In Refs.lTTl- 
I19l33l the master equation was obtained keeping the non- 
secular terms and this procedure produces the expected 
result. In a similar way, it can be verified that if we 
use the secular approximation'^ and perform the same 
calculation we obtain Tr5(£soc Q) 7^ 0. This result is in 
contradiction with the one based simply on a symmetry 
argument stating that since [Hj, Q] = 0, we do not have a 
dissipative contribution to the current. This observation 
is enough to state that the secular approximation is not 
feasible in the description of Co oper p air pumping as has 
already been pointed out ear lier.fi^Illl This result is much 
more general and can be discussed in completely abstract 
terms. 



B. Flux noise in a Cooper pair sluice 

If the Cooper pair sluice is subject to flux noise the dis- 
sipative contribution to the current is finite and, hence, 
requires further analysis. As an example we consider 
the case in which the fluctuations of the magnetic flux 
through the outer loop influence the total phase across 
the sluice ip. We will refer to this as phase bias noise. 

We consider small phase fluctuations in the vicinity of 



with 6ip <^ ipQ. The total Hamiltonian is Hj = Hj{ipo) + 
SHj with 



SHj = 



Jl sin 



(f-«)+.«.i..(f + ») 



Sip 

Y' 



(15) 



With the two-state approximation, we can write this 
in terms of the excess of Cooper pairs on the island 
{|0),|1)}. Using the formulas e*^ = |1)(0| and e"'^" = 
|0)(1|, we have 



dHj^{Sr\0){l\+dJ\l){Q\) 



(16) 



where SJ = (sin ^J++i cos ^ J-) and J± — Jl ± Jr- 



Equation ( 16 1 can be interpreted as the system- 



the static point ipo so that the total phase is 



2 



-6(p 



environment interaction Hamiltonian when 6(p are in- 
duced by the environment. Note that the operator 
acting on the system degrees of freedom is TJ/.s = 
(5J*|0)(1| + (5J|1)(0|. The charge operator on the is- 
land is Q = — 2e|l)(l| and, since [Hi_s,Q] 7^ 0, we 
must include the contribution of the dissipative current. 
As discussed above, this contribution can be calculated 
using the master equation approach. To this end, we 
can employ the general form of the master equation pre- 
sented in Ref. [T7]. The only difference resides in the 
matrix elements of the coupling operator (m|i7/^5|n), 
where \n) and |m) are the time-dependent eigenstates 
of the sluice Hamiltonian Hg, and in the spectral den- 
sity function of the environment. The latter is defined 
as Sip{uj) — dT{Sf{T)d(f{0))e'^'^'^ and it can be calcu- 
lated from the correlation function {5Lp{T)d(p{0)) of the 
phase fluctuations. 
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FIG. 2: (Color online) (a) The spectral density function of the phase bias noise as a function of the magnetic flux through 
the control SQUID (j!)/$o for frequency tJm = 1.7 X 10^° s-^ The inset is a zoom near the resonance point (j!)/$o = 1/2. The 
dashed lines show the region delimited by the minimum and maximum S,p{ujm, (j})- Its width is given by Lo/{-nL). (b) The 
ratio 5'^(cj,0)/5'^(tJ,O) for 0/$o = 0.5 (solid line) , 0/<I>o = 0.4991 (dashed line) and <^/$o = 0.5003 (dot-dashed line). The 
shadowed region denotes the frequency range spanned during a typical pumping cycle shown inFig.[l];b): 1.7 xlO^^ s'^ < a; < 
7.9 x 10^" s"^ The parameters of the artificial noise are 7? = 30 fi, i?s = 500 fi, Cs = 50 fF, Jc = 25 ^iA, L = M = 0.69 nH. 



IV. ENGINEERED ENVIRONMENT FOR 
PHASE BIAS NOISE 

To determine experimentally and discriminate the ef- 
fect of the environmental noise on the Cooper pair pump- 
ing, we should be able to control several properties of the 
environment. This is possible if, in addition to the natu- 
ral environment, the system is coupled to an engineered 
source of noise. 

A schematic description of the circuit used for the im- 
plementation of such a phase bias environment is pre- 
sented in Fig. [l] The main source of noise is a thermal 
resistor R with noise voltage 5V . This circuit is coupled 
to the system by mutual inductance M and thus it per- 
turbs the phase across the device. To modify the effects 
of the environment we introduce a control SQUID. By 
controlling the flux threading it, we can change the cur- 
rent noise in the circuit which is coupled to the system. 

The control SQUID can be represented by a RLC 
parallel circuit with resistance i?s, inductance Ls 
and capacitance Cs- The impedance of a sin- 



gle SQUID is Zrlc = 

iLsi<p)ujRs 



RsZcqZn 



■ r — , p — r ,,s \ , whcre Zr„ and Zn^ are the 
impedances associated to the inductance Ls and ca- 
pacitance Cs, respectively. The SQUID inductance 
Ls{<t>) — LqI cos(7r(/)/<I>o) can be controlled by chang- 
ing the normalized flux through the SQUID (p/^Q. The 
constant component depends on the maximum critical 
current of the SQUID Ic'- Lq = h/{2nelc)- Given 
Zrlc, can calculate the total impedance of the cir- 
cuit Ztot — Zrlc + Zr + Z^ where Zr and Zl are the 
impedances of the resistor and inductor, respectively. 

We consider a cold resistor with a cut-off frequency 
much higher than the typical transition energy of the 
system 51, i.e., fi ^ kgTR and ^ lOc where ujc, and Tr 
are the cut-off frequency and the temperature of the noise 
source, respectively. In this limit, the spectral density 
of the Johnson-Nyquist voltage noise across the resistor 
is Sv{uj) = 2hujR for > and Sv{<^) = for uj < 
0. We obtain the spectral density of the current noise 
5/(u;,0) = 5y(w)/|Ztotpas 



Si{uj,(f>) 



2RhLo 



R 



Ls(<l>)R%{l-CsLJ^Ls(<p)) I T 
R%{l-Csuj''Lsmf+Ls{4>yR% ^ 



(17) 



The flux noise spectrum is S'$(w, 0) = M'^Si{uj) and, 
hence, the phase noise spectrum reads 

S^{u:,(t))^ ^2 ■ (18) 



The behavior of S^{uj,(f>) as a function of the SQUID 
flux cj) and for a fixed frequency iVm is shown in Fig. [2ja) . 
The spectrum shows two resonances at (j)/^o = 1/2,3/2 
and it reflects the behavior of the control SQUID indue- 
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FIG. 3: (Color online) (a) The dynamic dissipative currents J^''*'='= (solid line) and the geometric dissipative currents l'^-'^^^^ 
(dashed line) normalized to the maximum geometric current I'j^ax during the pumping cycle. The maximum geometric current 
during one cycle is I^ax = 0.5 nA and the control SQUID of the environmental circuit is open with (;/)/<I>o ~ 1- (b) The 
pumped charge as a function of the flux through the control SQUID of the artificial environment circuit. In the inset we show 
the behavior near one of the resonance points. The dashed lines are determined from Fig. [2] (a) and they are located at the 
minimum — 1/2 and at <^/$o ~ 1/2 ± Lo/{nL) of the spectral density function. In the numerical simulations, the phase 

bias is (po — tv/2, the pumping frequency is / = 150 MHz, and the device parameters are jf'^ / Ec = 0.1, Ji^jjf^ = 0.03 (with 
i = L,i?), nf = 0.8 = 0.2, and Bc/ftB = 1 K \Ec I {I-kTi) = 21 GHz]. 



tance Ls{(j)). At these points the artificial environment 
is maximally decoupled. 

It is useful to analyze in detail the behavior of the 
spectrum near the resonance 4>/^o — 1/2. The spectral 
density function in Eq. ( |18[ ) for realistic circuit parame- 
ters used in the numerical simulations (see Fig. [2]) is well 
approximated by 



i?2 + L2^2 



1 - 



7rL(0/$o-l/2) 



V^max 



(19) 

This expression allows us to calculate the positions and 
the values of the maximum and the minimum of the 
spectral density function. The maximum is found at 
/$o = 1/2 + Lo/{7rL) and it is 5-'^^(w,„, <^n,ax) = 
The minimum can be obtained from Eq. 

(171 in the limit 0/$o ^ 1/2 and it is denoted by 
5^""(a;„,l/2). The width of the drop is given by the 
difference ^max/'&o — 1/2 = Lq/{'kL). It is determined 
by the ratio between the minimum inductance of the 
SQUID and the inductance of the circuit [see the inset 
in Fig. [2] (a)]. Using the circuit parameters in Fig. [2j 
we can estimate the decrease of the noise spectrum at 
the minimum: 5--(c<.„„ l/2)/5^-'^(c^„, 0,„ax)} « 10"* 
and S'™"(a;™,l/2)/S'^(a;„,0) « IQ-^. This means that 
we can reduce the strength of the system-environment 
interaction and effectively decouple the system from the 
artificial noise source. 

The noise influences the circuit at the frequency equal 
to the energy gap f2 in the system. In steered evolu- 
tion ri changes in time and so does the coupling fre- 



quency. To have an estimate of this effect, Fig. [2] (b) 
shows the change in the spectrum as the ratio between 
Sip{uj, (j))/Sip{uj, 0) for the range of frequencies f2 spanned 
in the evolution during the cycle in Fig. [T] (b) . Here con- 
trolling the SQUID reduces drastically the strength and 
the spectral density of the noise during the evolution. 

We note that the corresponding decoherence processes 
induced by the artificial noise are stronger than the ones 
naturally produced by the flux noise typically observed 
in the experiment (see Appendix [B|) . For this reason, we 
neglect the effects of the latter and only consider the 
decoherence induced by the engineered environment. 

In the following we also neglect the charge noise. An es- 
timate of the artificial flux noise strength (see Appendix 
[b| shows that, in average, it is stronger than the charge 
noise except that near the resonance points of Sip{uj,(j)). 
Thus, the effect of artificial flux noise near these regions 
would be partially hidden. To observe clearly the influ- 
ence of the flux noise along for all the values of (j), it 
is necessary to reduce the charge noise or to further in- 
crease the strength of flux noise coupling, i.e., increase 
the mutual inductance. 



V. NUMERICAL RESULTS 

Let us study the effect of the artificial noise on the 
pumped charge in the device shown in Fig. [ija). The 
SQUID energies and the gate voltage are modulated 
within the adiabatic time Tad = 1// where / is the pump- 
ing frequency. The system is initialized in the ground 
state and the dynamics is obtained by solving the master 
equation for the reduced density operator up to the first 
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order in the adiabatic parameterl^^Esl The pumping cy- 
cle is repeated until the system reaches the steady state 
solution ps[t + Tad) = ps{t)- We assume that the dis- 
sipative source R is at such a low temperature that the 
environment cannot excite the system. 

In Fig. |3][a), we show the time dependence of the dis- 
sipative dynamic and geometric currents induced by the 
phase bias noise. The curves correspond to the dissipa- 
tive currents flowing through the circuit when the system 
reaches the steady state solution. The dissipative cur- 
rents are about three orders of magnitude smaller than 
the total geometric current. The corresponding dissipa- 
tive pumped charge is almost zero. This is not surprising 
since we are analyzing the system in the steady state 
solution, and since the dissipative currents have no fa- 
vored direction, the corresponding charge averages out 
over many pumping cycles. The dissipative contributions 
are significantly greater during the first pumping cycles 
when the system is still far from the steady state. 

In Fig. [sjjb), we present the stationary pumped charge 
under the influence of the artificial flux noise. The tuning 
parameter is the flux of the control SQUID 0/$o which 
allows us to effectively change the system-environment 
coupling. As shown, the pumped charge is almost con- 
stant except near the resonance points (/)/$o = 1/2,3/2. 
The global behavior of the pumped charge and its drop at 
the resonance points can be physically explained keeping 
in mind that here the system-environment coupling is ef- 
fectively suppressed [Fig. [2ja)]. The stationary pumped 
charge is determined by two competing effects. The non- 
adiabatic transitions tend to excite the system while the 
relaxation induced by the environment tends to keep the 
system in the ground state. Since the ground and excited 
states pump in opposite directions, the first effect reduces 
the pumped charge while the second one stabilizes it near 
to the expected value of one Cooper pair per cycle .^^^1111211 
The stationary solution is reached after many pumping 
cycles when the two effects are in balance. Away from the 
resonance points, the spectrum is almost constant pro- 
ducing a constant pumped charge. At the resonances, 
the results in Fig. [Sj^b) reflect the fact that, since we 
are effectively decreasing the system-environment cou- 
pling, the non-adiabatic excitations become more im- 
portant. Correspondingly, the excited-state component 
is enhanced in the final steady state and we observe 
a smaller pumped charge. The maximum drop of the 
pumped charge is (Qmax - Qmin)/Qmax ~ 15% and the 
drop from the pumped charge away from the resonance 

(Qo - Qmi„)/Qo - 10% [Oo = Q(0 = 0)]. 

However, the details of the pumped current in the 
proximity of the resonances are different from the ones 
exhibited by the environment spectrum S^{uj,(j)). This 
reflects the fact that the pumped charge is an observ- 
able that can be influenced by many different effects. In 
particular, the pumped charge first increases near the 
resonance points before decreasing. It has clearly s non- 
trivial behavior [inset in Fig. |3](b)]. One possible ex- 
planation of this behavior can be the influence of the 



Landau-Zener-Stiickelberg (LZS) interference.!^ During 
the pumping cycle in Fig. [ijb), the system crosses two 
avoided crossing where Ug ~ 1/2. Near these points 
Landau-Zener transitions occur and the phase difference 
between the excited and the ground state accumulated 
in the intermediate region leads to interference effects. 
The overall result is to change the ground-state popula- 
tion and thus the pumped charge. In the present model, 
the phase bias is kept constant but the artificial noise 
can effectively induce a change in the accumulated phase 
difference and hence in the LZS interference. If the sys- 
tem is initially in the ground state, the probability to 
excite it due to LZS interferometry after one cycle is 
Ve = Plz{1 - Plz) cos2(a + (^/2) where Plz is the 
Landau-Zener transition probability and a depends on 
the energy gap and on the pumping frequency.^^ If 5if is 
the perturbation on the constant phase bias ~ 7r/2, 
we have (p/2 = ir/A + 5(p. The average of the excitation 
probability during a loop is 

'Pe ~ Vl + \plz{1 ~ Plz) sin(2a)( V), (20) 

where V^^ is the excitation probability without noise, 
((5<y9^) is the average of the square phase fluctuation and 
we have assumed (5<p) = 0. We calculated the contribu- 
tion due to the noise as {Sip^) = S^{uj,(j))dui and, 
as a function of the control flux </>, and observe that it 
behaves as S^p{uj,(f)) as a function of cj). The probability 
to remain in the ground state is Vg ~ 1 ^ Ve and its 
dependence on the control SQUID flux (p is inverted with 
respect to S{uj, 0). This means that when the strength of 
the noise spectrum increases because of the LZS effect, 
the ground state is less populated and the pumped charge 
decreases. The final behavior of the pumped charge as 
a function of is determined by the noise spectrum, the 
LZS interference and their relative strength integrated 
over the pumping cycles. 

Note that for the numerical simulations we have 
used the master equation which includes only the 
first order con tribution in the perturbative adiabatic 
parameter .i^^ESl Xo obtain a more accurate value of the 
pumped charge a different master equa tion which in- 
cludes high-order corrections can be used.li^lMl 



VI. CONCLUSIONS 

In summary, we have presented a study of the fiux 
noise on the Cooper pair pumping process. Firstly, we 
showed that in the presence of flux noise the current op- 
erator must be modified in order to take into account the 
current induced by the environment. This dissipative 
current is related to the symmetry of the system and can 
be calculated from the master equation for the reduced 
density matrix. 

Secondly, we then analyzed the effects of the phase 
bias noise on the Cooper pair pumping. In the model 
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studied, the noise is produced by an artificial environ- 
ment coupled to the system by mutual inductance. The 
advantage of this scheme is that, introducing a control 
SQUID in the environmental circuit, we have access to 
the system-environment coupling and can essentially de- 
couple the system from the environment. The pumped 
charge obtained by solving the master equation shows 
clearly the features induced by the noise. 

The system presented here operates in the charge 
regime and it is then primarily sensitive to charge noise. 
For this reason, it can be challenging to experimentally 
measure the predicted effect because it can be partially 
hidden by the charge noise. However, there are several 
modifications which can increase the experimental ac- 
cessibility: reducing the effect of charge noise increasing 
the mutual inductance coupling with a different design 
of the artificial noise circuit. For example, using mul- 
tilayer lithographic techniques, it is possible to increase 
the coupling between the system and the environment. 

This work is a step towards understanding the effect 
of different types of noises on steered superconducting 
devices. Such understanding is critical for practical ap- 
plications such as metrology and quantum information 
processing. A deeper knowledge of the effect of the en- 
vironment can open a way to the design and implemen- 
tation of robust devices or, in the case of quantum infor- 
mation, error correction techniques. 

In this direction, the pioneering experimental works 
of implementing engineered environments have been car- 
ried out using trapped ions.'^* One of the most strik- 
ing results of these experiments was the proof that it is 
possible to store quantum information in states which 
are robust against decoherence.l^ Still, apart from a few 
theoretical proposals,"^ similar experiments are missing 
in condensed matter systems. 

Along these lines, the next step is to analyze the effect 
of similar flux noise sources on superconducting devices 
in different regimes. For example, in the transmon regime 
the system is insensitive to charge noise and the flux noise 
becomes the dominant source of noisc:^'^ Thus, using the 
discussed engineered environment, we should be able to 
study and measure the environmental effect in controlled 
situations. 
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Appendix A: Absence of dissipative current with 
charge noise environment 



By using Eq. ([s]) we can check by a direct calcula- 
tion that, when the system is subjected to charge noise, 
there is no dissipative current contribution. When the 
system is adiabatically steered along a cyclic loop, the 
master equation is conveniently written in the time de- 
pendent eigenbasi s of J?g. There are several equiva- 
lent approacheJ^^^^^'^ exploiting the two-level approx- 
imation in the weak coupling limit. Since the present 
calculation can be done in an abstract way without ex- 
plicitly taking into consideration the form of Hs, it is 
convenient to use the ma ster e quation in the superadia- 
batic basis representatiorP^l^s] which is easier to manip- 
ulate. 

Let Di be the operator which diagonalizes Hs in a 
fixed basis. Then the dynamics of the closed system are 
effectively governed by Hs.i — d\HsDi + hwi where 
wi ~ -~iD\Di. Defining D2 as the diagonalizing oper- 
ator for Hs,i allows us to define a new effective Hamil- 

tonian as Hs,2 = D\HsD2 + tvw2 where W2 = ^iD\D2- 
If the evolution is sufficiently slow, that is, \'W2\ ^ \wi\ 
we may approximate Hs,2 ~ D\HsD2 yielding an effec- 
tively nonsteered picture. Denoting the superadiabatic 
basis as {\g) = ^li^alO), |e) = DM^)), where {|0), |1)} 
defines the fixed basis, yields a master equation in the 
Schrodinger picture as 



= -(r^e + '^eg)Pgg + SRclfoPge} + ^eg, (Al) 



and 



= iVtpge - (f + + T-)Pgg 

(A2) 



where Prs = {t\ps\s) with r, s e {<?, e}, and il = [E^. — 
Eg)/h. If Hi = Z (S) X where Z and X are the operators 
acting on the system and the environment, respectively. 
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the rates in Eqs. ([Al|) and (|A2|) are defined as 



MZ\g)\ 



\{e\Z\gW 



r — 

-L eg — 



~ {e\Z\g}i{g\Z\g) - {e\Z\e)) 

io - i'(Oj, 



{g\Z\e){{e\Z\e)-{g\Z\g)) 



\{e\Z\e)\^ , \{9\Zh 



s{±n), 



2h^ 



+ 



{9\Z\g){e\Z\e) 



2h^ 

{9\ZW 
{9\Z\e)^ 



(A3) 
S{0), 



s{n), 



2h^ 



The spectral density is denoted by SUu) = 
J^^dT{X{T)X{0))e"^^. Note that Eqs. (aiI and 



( |A2[ ) still contain the non-secular terms. 

We can write the superadiabatic states in a general 
form as 



\g) = cos6l|0) +sin6le*'^|l), 
|e) = sin6i|0) -cos6le-*'^|l) 



(A4) 



where 9 and 7 are time dependent functions. We con- 
sider a generic diagonal noise operator in the charge ba- 
sis: Z = a|l)(l| + b\0){0\ for a,b e R. Expressing the 
charge operator of the island as Q = — 2e|l)(l| and the 
noise operator Z in the time-dependent basis (A4), we 
have 



and 



Q = -2e[sm'^e\g){g\+cos^0\e){e\ 
-cos0sin0(|e)(,g| + |g)(e|)], 

Z = {acos^0 + bsm'^e)\g){g\ 
+{bcos'^e + asin^ 0)\e){e\ 
{a-b)cos0siii9i\e){g\ + \g){e\). 



(A6) 



At the same time the Lij terms in Eq. ([7| can be read 
directly from the dissipative part of the master equation 



^99 = -(r 



^eg)Pgg + ^e{ToPge} + Teg 

2 +^ 



f-)Pgg 



eg 
2 



r/3)Pe 



(A7) 



Thus we have all the elements to calculate the dissipative 
current 

Trs(£ Q) = {Qgg - Q,e)Cgg + 2^eiCgeQeg). (A8) 

By explicitly writing the transition rate in Eq. ( [A| for 
the noise operator Z and inserting all the contributions 



in Eq. (A8), we verify that Trs(£ Q) = 0. These results 



are valid for any perturbation operator diagonal in the 
{|0), |1)} basis and thus for the charge noise discussed in 
Sec. |III A| The same result can be obtained using the 
master equation in the adiabatic basis^^ 

It is interesting that the condition Tr5(£ Q) — de- 
pends critically on the form of the master equation. As 
an example, we consider the same problem but we per- 
form the usual secular approximation.!^ The correctness 
of the secular approximation in th e analysis of the charge 
pumping has been qucstionecP^E^ since it has been shown 
that it can lead to unphysical results such as charge non- 
conservation. The master equation after the secular ap- 
proximation reads 



dp] 



sec 

'99 



and 



,scc 



dt 



dt 



= -(r 



^ eg 

2 



^eg)Pgg + Teg) 



|_ge 
2 



^ <p I Pge ■ 



(A9) 



(AlO) 



In this case, it can be shown that a calculation similar 
to the one above gives non-vanishing dissipative current 
since Tr5(£sec Q) 7^ 0. This result is in contradiction 
to the one based on the symmetry argument. Similar 
results can be obtained in a more formal framework as 
discussed in Ref. [501 



Appendix B: Artificial and natural decoherence rates 

In addition to the artificial environment discussed in 
Sec. |IV[ we should also take into account the effect 
of the natural flux noise. It is sufficient to consider 
the generic decoherece rate r(0) = \{i\SH\j) / hl"^ S^{uj , 4>) 
(A5) where Sip{uj,(j)) is the spectral density function of the 
environment. The matrix element |(«|(5J?|j)|2 is associ- 
ated to the relaxation process if i ^ j and the dephasing 
if i = j. The order of magnitude of \{i\6H\j) /h\'^ for 
the pum ping cycle in Fig. [ijb) can be estimated using 
Eq. (16) as {Jm/K)^- Combining this with the values of 
we have an estimate of the decoherence times 



(both dephasing and relaxation), r((/)) = 1/T{(j)). 

During the pumping cycle the system energy gap 
changes in tim d^^ l ^^ and the minimum energy gap is 
reached near the degeneracy points where Ug = 1/2. For 
the majority of the evolution time, the system frequency 
is close to its maximum ujm = 7.9 x 10^° s~^. We can 
use this reference frequency to estimate the decoherence 
rates. The minimum and maximum values of the spec- 
tral density function at 00 m are S"™™ = 2.8 x 10~^^ s and 



5" 



3 X 10 s. The estimated relaxation and de- 



pahsing times with the parameters used in the numerical 
simulation (see the caption of Fig. [3]) are Tmin ~ 1.9 ns 
and Tmax ^2 us. 

For the charge noise, the decay and dephasing times to 
the can be estimated as in Ref. [T7] and they are of order 
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10 ns. Thus, the artificial flux noise is the dominant 
source of noise except that near the resonance points of 

The effect of the natural environment can be estimated 
from the experimental results of the decoherence time 
for devices affected by flusP^"^. We focus on the ef- 
fect of ubiquitous 1// noise. In this case, the measure- 
ments of the decoherence induced by such a noise sug- 
gest that the decay has Gaussian and not exponential 
shap^^2143| Note that, strictly speaking, this implies that 
it cannot be described by a linear master equation. How- 
ever, a phenomenological approach has been appliecP^. 
The effective dephasing rate is F = \/ ^n2\drt/d(p\ 
where the spectral density function for the 1// noise 
is S^{uj) — A^p/\ll)\ and Q is the frequency of the sys- 



tem. The maximum of \d^/d(p\ is obtained when the 
system is at the degeneracy point, i.e., ng = 1/2 and only 
the Josephson contribution is present Hamiltonian [see 
Eq. In this case, n = J% + Jl + ^JrJl cos ip/h 

and we can approximate Jr ~ Jj^.j ^ « J^. From 
the measured amplitude of the magnetic flux spectrum 
A$ = (1.7 /z$o)^, we can calculate the amplitude of the 
phase spectrum = Att^A^/^q. Thus, we obtain the 
dephasing time due to low-frequency noise: r « 28 /iS. 
This dephasing time is very long compared to the ar- 
tificial decoherence time, and hence the natural low- 
frequency noise does not change the results of the pre- 
sented numerical simulations. 
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